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Abstract. In this talk we discuss the possibility of constructing a fluctuation theory 
for structural glasses in the non-equilibrium aging state. After reviewing well known 
results in a toy model we discuss some of the key assumptions which support the validity 
of this theory, in particular the role of the configurational entropy and its relation to 
the effective temperature. Recent numerical results for mean-field finite-size glasses 
agree with this scenario. 



I THE GLASS STATE 

vb. 

A theoretical understanding of the physical mechanisms behind the glass transi- 
tion remains an open problem [1]. What is commonly understood as a glass (for 
instance, window glass) is a metastable phase, with free energy higher than that 
of the crystal obtained by the continuation of the liquid line below the melting 
transition temperature T M . Such a glassy state may be experimentally achieved by 
cooling the liquid fast enough. Glasses share physical properties common to liquids 

_L ', and solids making difficult to decide in which phase (if any) their are. On the one 

hand, atom positions are randomly located in the glass much alike a liquid and 
apparently there is no long-range order. On the other hand, under compression 
glasses behave like a solid showing a very low mobility. Glasses constitute a state 
of matter in between liquids and solids, not completely classificable as any of them, 

rS ' which show the following generic features: 

• The viscosity anomaly. Glasses show a very rapid increase of the viscosity 77 
in a relatively narrow range of temperatures. For instance, in window glass 
the viscosity increases over nearly 20 orders of magnitude by changing the 
temperature in a range within a 10% of the value of the melting transition 
temperature. This increase is often well fitted by the Vogel-Tamman-Fulcher 
law, i] ~ exp(A/(T — T )) where T and A are free parameters. Fragile glasses 



are those where T is finite while for strong glasses T is small (compared to 
Tm) and the viscosity displays Arrhenius behavior. By convention the glass 
transition temperature T g is defined such that r)(T g ) = 10 13 Poise which for 
many liquids corresponds to a relaxation time of several minutes . 

• Two-step relaxation. In the undercooled regime glasses show relaxation func- 
tions with a characteristic two-step form. Intermediate scattering functions 
show a first decay to a plateau (/3-process) followed by a secondary and slower 
relaxation (a-process) which defines the longest activated time. Although both 
process are activated their characteristic times turn out to be well separated 
below the mode-coupling transition temperature. 

• Aging. Suppose a glass is quenched to a temperature Tf below the glass 
transition temperature T g . If one were able to measure the viscosity as a 
function of time one would observe that it grows with time t according to the 
approximate law: 77 ~ 770 + at where r]o is the initial value reached soon after 
the quenching and a is a temperature dependent parameter a ~ exp(— B /Tf) 
where B is an activation barrier. This growth can be extremely slow since 
a can be very small depending on the value of the ratio B /Tf. The growth 
of the viscosity manifests itself as a dependence of intermediate scattering 
functions on the time of the measurement (usually called waiting time t w ) 
where the longest decorrelation time depends on the value of the viscosity rj 
at t w . The dependence of the response of the system in the presence of an 
external perturbation, which softens as time goes by, is commonly referred as 
aging. 

II A SIMPLE SOLVABLE MODEL OF A GLASS 

Before describing the key assumptions to understand glasses it is convenient to 
consider an instructive example. One of the simplest solvable models which shows 
key features of structural glasses is the oscillator model with parallel dynamics 
introduced by one of us in [2]. This model contains some essential features of non- 
equilibrium thermodynamics applied to glassy systems. The model is defined by a 
set of non-interacting N harmonic oscillators with energy: 

K N 
E = — Y^ tf, - 00 < Xi < 00 (1) 

^ i=i 

where K is a coupling constant. An effective interaction between oscillators is 
introduced through a parallel Monte Carlo dynamics characterized by small jumps 
Xi — > x\ = Xi + -?=; 1 < % < N where the variables r,i are extracted from a Gaussian 

distribution P(r) = , 1 exp(— r 2 /(2A 2 )). At each Monte Carlo step all oscillators 
are updated following the previous rule and the move is accepted according to 
the Metropolis algorithm with probability W(AE) = min[l,exp(— /3A-E 1 )] where 



AE = E({x'}) — E({x}). The smallness of the jumps in the variables Xi is required 
in order to give a finite change in the energy so the acceptance does not vanish in 
the N — ► oo limit. 

This model (as well as some modifications proposed afterwards [3]) has been 
extensively studied and shows the following scenario. Because of the finiteness of the 
moves the ground state cannot be reached in a finite amount of time. Despite the 
absence interactions in the Hamiltonian, the Monte Carlo dynamics induces entropy 
barriers corresponding to the flat directions in energy space that the system hardly 
finds when the acceptance is low. A simple calculation at finite temperature shows 
activated behavior for the relaxation time r re i ax ~ exp(^-) despite the absence 
of energy barriers in the model (the potential is a single well in N dimensions) 
making the dynamics of this model quite reminiscent (but simpler) of that of the 
Backgammon model [4]. 

The interesting dynamics in this model is found when studying the relaxation 
after quenching to zero temperature. In what follows we summarize the main 
findings for this case: [2]: 

• Slow decay of the energy: The evolution equation for the energy is Markovian. 
This simplicity allows for an asymptotic large-time expansion. The energy 
asymptotically decays logarithmically E(t) ~ l/log(t) and the acceptance 
ratio decays faster A(t) ~ l/(t log(t)). 

• The effective temperature. The fluctuation-dissipation (FDT) ratio [5] can be 
exactly computed and depends only on the lowest time. At zero temperature 
one gets in the large s limit, 

dC(t,s)/ds 2E(s) 

Tcff(s) - G(t, 3 ) ~ir (2) 

i.e., equipartitioning is obeyed off-equilibrium allowing to define an effective 
temperature in terms of the FDT ratio. 

• The role of the configurational entropy. The effective temperature previously 
obtained allows for a definition of a time dependent configurational entropy. 
At time t the number of configurations explored by the system is given by the 
surface of an iV-dimensional sphere of radius R = E? where E is the energy 
at time t. The configurational entropy is then given by Q ~ R N ~ X = E~z~ 
leading to an extensive configurational entropy S c = log(fi) = y log(E) which 
satisfies, using (2), the canonical thermodynamic relations, 

dE 2E T eff l ' 

• Heat-flow driven by the effective temperature. Following [6] we can obtain the 
resulting heat-current by coupling two identical harmonic oscillator systems 
with a term ex^xji and measuring the energy flow to order e 2 . One gets [7], 



J(t) = hm(^M - ^M) (4) 

t^ s at as 

with Q(t,s) = jjYliLi x i(t)Ui( s )- If one °f the two systems acts as a ther- 
mometer with characteristic frequency cv >> 1/t then the measured temper- 
ature coincides with the effective temperature and the Onsager relation is 
satisfied, 

J = L QQ V{^) (5) 

with Lqq ~ l/(£ log(t)) and A ~ log(t)/t for the thermal conductivity [7]. 
This example generalizes the Fourier law (and by extension, all linear relations 
between fluxes and forces initially derived in the vicinity equilibrium [8]) to 
the case where fluctuations are present around the off-equilibrium aging state. 

Ill THE CANONICAL THERMODYNAMIC 
SCENARIO FOR GLASSES 

In the previous example we have seen how equipartitioning in the off-equilibrium 
state is satisfied when the acceptance rate is small. This can be easily understood: 
at zero temperature only configurations with the same or lower energy are accepted. 
Because the size of the moves A is finite the system samples all the configurations 
of equal or lower energy with the same probability. Being their number a mono- 
tonically increasing function of the energy and because the energy cannot increase 
at T = we conclude that the relevant fluctuations are those which explore the 
constant energy surface. After quenching to T = and for long times, the prob- 
ability to explore a different state with energy E < E* = E(t) is uniform for all 
states with a given energy and given by, 

P ~ n(E) = exp(NS c (E)) ~ exp l — log(E) 

exp (^ \og(E*) + ^(E- £*)) oc exp((3 cS (E - E*)) &(E* - E) . (6) 

Fluctuations are then described in terms of an effective temperature given by the 
configurational entropy as described in eq.(3). In the last years Th. M. Nieuwen- 
huizen, inspired by results for this oscillator model as well as for the p-spin spherical 
spin glass has proposed that a similar scenario could apply for generic glassy sys- 
tems [9]. Although it is not clear in what conditions these results are valid for 
realistic systems the idea is suggestive enough to merit a more detailed investiga- 
tion. A more elaborate picture along the same line has been presented by Franz and 
Virasoro who, using well known results from spin glasses and assuming a separa- 
tion of timescales, have proposed that a modified version of the Onsager hypothesis 



should put the configurational entropy into the game [10]. Although their analysis 
is based on some mean-field models it can be probably extended to realistic systems 
under some general mild assumptions. 

The idea goes as follows. Imagine a glass that is obtained by quenching the liquid 
to a temperature Tf below T g . As already discussed, one of the main experimen- 
tal facts regarding relaxational processes (for which mode-coupling theory gives a 
fairly good description in a large range of temperatures) is the emergence of well 
separated timescales (the a and f3 processes). Dynamics can then be viewed as if 
the system were jumping among different metastable states or basins, each basin 
defined as the set of configurations explored by the system during a time t* of the 
order of the relaxation time itself. The definition of the precise value of t* does 
not matter as soon as the timescales separation becomes very strong, a condition 
which emphasizes the robustness of the present scenario. Fluctuations around a 
given configuration have two contributions: 1) Thermal fluctuations within a basin, 
i.e., typical fluctuations occurring on a timescale less than t* and driven by the tem- 
perature Tf and 2) Activated jumps between different basins typically occurring 
for times larger than t* . If basins with a given free energy are more or less similar 
(they are not too big or too small) then their typical lifetime is of order t* implying 
that the probability to visit them at subsequent times t > t* is simply proportional 
to their number exp(S c (F)) where S C (F) defines the configurational entropy, 

P(F) oc fi(F) = exp(S c (F)) = exp(S c (F*) + (?^1) F=F .(F - F*) 

,.x 2 \ (,dSc{F), 



+0(F - F*Y) oc exp((^^) F=F »(F - F*)) (7) 

where F* is the threshold value of the free energy, i.e., the free energy of the basins 
at time t. Note that the present threshold we define here is conceptually different 
from that defined in mean-field p-spin models since now F* is a time dependent 
quantity. From (7) we can infer that the probability distribution for basin-to- 
basin fluctuations is given by the configurational entropy S C (F) which satisfies the 
relation, 

J_ = ( dScjF) , 

T cS v OF )f=f* k > 

and defines an effective temperature for interbasin fluctuations. The essentials of 
the physics contained in eq.(7) is the same as that contained in eq.(6) for the oscil- 
lator model with the difference that in the oscillator model fluctuations above the 
threshold E* are forbidden because Tf — while now fluctuations are dominated 
by jumps to basins with free energies around the threshold F* and differing by a 
finite amount. We must stress that eq.(7) gives a good description of the fluctua- 
tions when Tf is much smaller than T e g, otherwise thermal activation driven by Tf 
will modify the simple entropic contribution entering in (7). Biroli and Kurchan 
[11] have recently addressed the notion of configurational entropy in the framework 



of the Gaveau-Schulmann formalism showing that S c is a time-dependent quantity. 
This is in complete agreement with the present scenario: basins with free energies 
F above the time-dependent threshold F* are completely washed out as time goes 
leading to the appearance of a time-dependent cutoff S C (F) = S C (F) F < F* and 
S C (F) = — oo F > F* with S C (F) a monotonous increasing function of F. This 
keeps the full meaning of (8). 

The conclusion of this discussion is that, as soon as the formula (7) gives a faithful 
description of basin-to-basin fluctuations occurring beyond a timescale t* (which 
changes itself as a function of time), then it is possible to write down a dynamical 
free energy which controls both fast and slow processes: 



dyn 



Fbasin — T e ff S C (F) . (9) 



Keeping in mind that both terms contribute to physically different processes occur- 
ring in different (well separated) timescales this corresponds to the formula obtained 
in the framework of p-spin models. We must stress that in models where free energy 
basins are not uniformly sampled the simple results (8,9) cannot hold. Examples 
of such a class of systems are generic coarsening models into a finite number of 
absorbing states [12]. 

This formula has been recently checked in finite-size mean-field models for glasses 
where activated processes occur with finite probability [13,14]. One can evidenciate 
the presence of an activated regime different from that observed in the N — ► oo 
where ergodicity breaks below the threshold. To verify the above scenario is neces- 
sary to estimate the configurational entropy. A good recipe to do that is to partition 
the phase space into basins by counting the number of energy minima. These are 
also called inherent structures and were proposed by Stillinger and Weber as a 
powerful tool to investigate the landscape energy surface of liquids [15]. The above 
scenario has been verified in the ROM model [16] which is a faithful microscopic 
realization of the random energy model introduced by Derrida [17]. In this model, 
Fb as i n — Ebasin because the intra-basin entropy is very small (qEA ~ 0.96 for the 
metastable states) and the effective temperature can then be obtained using eq. 
(8) with S C (F) ~ S C (E). The results for the configurational entropy and the ef- 
fective temperature (obtained through FDT plots) are shown in figure (1), panels 
(a) and (b) respectively. The agreement between equation (8) and the numerics 
is excellent. Recent results [18] in the case of Lennard- Jones glasses, where the 
contribution Fbasin turns out to be important, also confirm these results. 

The simplicity of the hypothesis presented here suggests that some generic and 
simple principles are behind the statistics of fluctuations around the aging state for 
structural glasses. Surely we will see fast developments in the forthcoming years 
which will show to what extent these hypothesis are correct. 

Acknowledgements We are grateful to A. Garriga for a careful reading of the 
manuscript. 
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FIGURE 1. (a) Configurational entropy as a function of energy. The data are temperatures 
T = 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. The line is the quadratic best-fit. From Rcf. [13]. (b) FDT 
plot for the ROM model (Integrated response function as a function of IS correlation function) . 
The dash line has slope (3{ = ^- = 5.0, while the full lines is the prediction S C (E) from Ref. [14]: 
T cf f(2 n ) ~ 0.694, T cff (2 16 ) ~ 0.634 and T cff (2 19 ) ~ 0.608. The dot-dashed line is (3 cS for t w = 2 11 
drawn for comparison. 
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